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' ' Abstract 

^: 

The confinement mechanism earlier proposed by author is applied to describe the 
(possible) ground state of toponium r\t- For this aim the nonperturbative consis- 
tent approach is elaborated in both relativistic and nonrelativistic cases. The study 
entails estimates for parameters of the confining SU(3)-gluonic field in the above 
quarkonium, those estimates being also consistent with possible width of decay 
rjt — > 27. The corresponding estimates of the gluon concentrations, electric and 
magnetic colour field strengths are also adduced for the mentioned field at the 
t^T* ' scales of toponium. 
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In Refs. [1-3] for the Dirac- Yang-Mills system derived from QCD-Lagrangian 
there was found and explored an unique family of compatible nonperturbative 
solutions which could pretend to decsribing confinement of two quarks. Ap- 
plications of the family to description of both the heavy quarkonia spectra [4] 
and a number of properties of pions, kaons and //-meson [5-7] showed that the 
confinement mechanism is qualitatively the same for both light mesons and 
heavy quarkonia. At the moment it can be decribed in the following form. 

Two main physical reasons for linear confinement in the mechanism under dis- 
cussion are the following ones. The first one is that gluon exchange between 
quarks is realized with the propagator different from the photon-like one and 
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existence and form of such a propagator is direct consequence of the unique 
confining nonperturbative solutions of the Yang-Mills equations [2,3]. The sec- 
ond reason is that, owing to the structure of mentioned propagator, quarks 
mainly emit and interchange with the soft gluons so the gluon condensate (a 
classical gluon field) between quarks basically consists of soft gluons (for more 
details see Refs. [2,3]) but, because of that any gluon also emits gluons (still 
softer), the corresponding gluon concentrations rapidly become huge and form 
the linear confining magnetic colour field of enormous strengths which leads 
to confinement of quarks. This is by virtue of the fact that just magnetic part 
of the mentioned propagator is responsible for larger portion of gluon con- 
centrations at large distances since the magnetic part has stronger infrared 
singularities than the electric one. Under the circumstances physically nonlin- 
earity of the Yang-Mills equations effectively vanishes so the latter possess the 
unique nonperturbative confining solutions of the Abelian-like form (with the 
values in Cartan subalgebra of SU(3)-Lie algebra) [2,3] that describe the gluon 
condensate under consideration. Moreover, since the overwhelming majority 
of gluons is soft they cannot leave hadron (meson) until some gluon obtains 
additional energy (due to an external reason) to rush out. So we deal with 
confinement of gluons as well. 

The approach under discussion equips us with the explicit wave functions for 
every two quarks (meson or quarkonium). The wave functions are parametrized 
by a set of real constants dj,bj,Bj describing the mentioned nonperturbative 
confining gluon field (the gluon condensate) and they are nonperturbative mod- 
ulo square integrable solutions of the Dirac equation in the above confining 
SU(3)-field and also depend on /x , the reduced mass of the current masses of 
quarks forming meson. It is clear that under the given approach just constants 
a,j, bj, Bj, Ho determine all properties of any meson (quarkonium), i. e., the ap- 
proach directly appeals to quark and gluonic degrees of freedom as should be 
according to the first principles of QCD. Also it is clear that the mentioned 
constants should be extracted from experimental data. 

Such a program has been to a certain extent advanced in Refs. [4-7]. The 
aim of the present paper is to return to heavy quarkonia physics to obtain 
estimates for aj, bj, Bj in the possible ground state of the heaviest quarkonium 
- toponium (in what follows we denote it as r) t ), where so far little is still known 
about experimental spectroscopy of the system. 

Under the situation a certain motivation of studying toponium is adduced 
in Section 2. Section 3 contains main relations underlying both relativistic 
and nonrelativistic descriptions of any quarkonia in our approach. Section 
4 is devoted to computing electric form factor, the root-mean-square radius 
< r > and magnetic moment of the quarkonium under consideration in an 
explicit analytic form. Section 5 gives an independent estimate for < r > which 
is used in Section 6 for obtaining estimates for parameters of the confining 
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SU(3)-gluonic field for the toponium ground state rjt in both relativistic and 
nonrelativistic cases. Further in Section 7 we show that estimates of Section 
6 can also be consistent with possible width of two-photon decay r\ t — > 27. 
Section 8 employs the obtained parameters of SU(3)-gluonic field to get the 
corresponding estimates for such characteristics of the mentioned field as gluon 
concentrations, electric and magnetic colour field strengths at the scales of i] t 
while Section 9 is devoted to discussion and concluding remarks. 

At last Appendices A and B contain the detailed description of main building 
blocks for meson wave functions in the approach under discussion, respectively: 
eigenspinors of the Euclidean Dirac operator on two-sphere § 2 and radial parts 
for the modulo square integrable solutions of Dirac equation in the confining 
SU(3)-Yang-Mills field while Appendix C supplements Appendices A and B 
with proof of the fact that the so-called nonrelativistic confining potentials do 
not obey the Maxwell or SU (3)- Yang-Mills equations. 

Further we shall deal with the metric of the flat Minkowski spacetime M that 
we write down (using the ordinary set of local spherical coordinates r, 1?, ip for 
the spatial part) in the form 

ds 2 = g^da? ® dx u = dt 2 - dr 2 - r 2 (dd 2 + sin 2 tidtp 2 ) . (1) 

Besides, we have \8\ = | det(^ AU/ )| = (r 2 sin^) 2 and < r < 00, < < 7r, 
< if < 2tt. 

Throughout the paper we employ the Heaviside-Lorentz system of units with 
H — c — 1, unless explicitly stated otherwise, so the gauge coupling constant g 
and the strong coupling constant a s are connected by relation g 2 /(47i) = a s . 
In what follows we shall denote L 2 (F) the set of the modulo square inte- 
grable complex functions on any manifold F furnished with an integration 
measure, then L%(F) will be the n-fold direct product of L 2 (F) endowed with 
the obvious scalar product while t and * stand, respectively, for Hermitian 
and complex conjugation. Our choice of Dirac 7-matrices conforms to the so- 
called standard representation and is the same as in Ref. [5]. At last <g> means 
tensorial product of matrices and /„ is the unit n x n matrix so that, e.g., we 
have 

/V \ 

7 3 ® y = 7^ 

V 7^ / 

for any Dirac 7-matrix 7 M and so forth. 

When calculating we apply the relations 1 GeV -1 w 0.1973269679 fm, 1 s _1 w 
0.658211915 x 10~ 24 GeV , 1 V/m w 0.2309956375 x 10~ 23 GeV 2 , 1 T = 
4tt x 10- 7 H/m x 1 A/m w 0.6925075988 x 10^ 15 GeV 2 . 

Finally, for the necessary estimates we shall employ the T 00 -component (volu- 
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metric energy density ) of the energy-momentum tensor for a SU(3)-Yang-Mills 
field which should be written in the chosen system of units in the form 

T,u = ~F; a <f# + -F^ F^g^g^ . (2) 

2 Motivation 

Though theoretical toponium physics has been developing already during a 
long time (see early reviews of Refs. [8] and references therein) so far little is 
still known about experimental spectroscopy of the system, as we have men- 
tioned in Section 1. According to standard model (SM) with three generations 
(see, e.g., Ref. [15]) main properties of top-quark are determined by summand 
in SM-Lagrangian of the form —g/(2y/2) [£7^(1 + 75) V#, where the gauge 

coupling constant g is connected with the Fermi coupling constant Gp as 
g 2 \/2 = 8G F rriw (m w xs 80.403 Gev is mass of lU-boson) while V tb is the 
corresponding element of the Cabibbo-Kobayashi-Maskawa mixing matrix. 
In second order in g this yields the decay width T(t — > W b) ~ GV|Vtf,| 2 mj? so 
that T proves to be between (1.0 - 1.6) GeV depending on top-quark mass m t 
[15]. With its respectively short lifetime of order 0.5 x 10~ 24 s, the top quark is 
expected to decay before top-flavored hadrons or tt-quarkonium bound states 
can form although as far back as in Ref. [9] possibilities for formation of to- 
ponium at m t ~ 170 GeV were discussed (see also Refs. [10] and references 
therein) . 

The above estimates, however, suppose \V t b\ to be of order 1. But it is known 
that mixings and the number of fermion generations are not fixed by SM. Un- 
der the circumstances possible existence of extra SM families may sufficiently 
decrease \ V t b\ so that toponium can be formed [11]. Under this situation, one of 
the main signals for detecting toponium should be decay r\ t — > 27 and enough 
number of those decays might be observable even at LHC and its further 
upgrades and, of course, at possible future collider VLHC (for more details 
see Ref. [11]). For the sake of justice we should, however, note that present 
measurements still indicate that \V t b\ >0.78 [15]. 

Also within the SM based on the two Higgs doublets formation and observation 
of heavy quarkonia (including toponium) might be lightened [12]. Thus, on the 
whole, the question of existence and observability for toponium remains open 
and, probably, possible discovery of toponium at future colliders might provide 
a connection to new or unexpected physics. 

Specifically, therefore, there is an certain interest of trying to explore the topo- 
nium already now from the point of view of the above confinement mechanism 
which directly appeals to quark and gluonic degrees of freedom as should be 
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according to the first principles of QCD. Under the circumstances we shall 
use some previous theoretical estimates for toponium obtained from other 
considerations (see, e. g., Refs. [13,14] and references therein). 



3 Main relations 



As was mentioned above, our considerations shall be based on the unique fam- 
ily of compatible nonperturbative solutions for the Dirac- Yang-Mills system 
(derived from QCD-Lagrangian) studied in details in Refs. [1-3]. Referring for 
more details to those references, let us briefly describe and specify only the 
relations necessary to us in the present paper. 

One part of the mentioned family is presented by the unique nonperturbative 
confining solution of the Yang-Mills equations for A = A^dx^ = A^Xadx^ 1 (X a 
are the known Gell-Mann matrices, /i — t,r,~&,(p, a — 1, ...,8) and looks as 
follows 



A lv = Al + -^=A% = hr + B 1 , A 2(p = -A% + -j=A% = b 2 r + B 2 , 

A 3v = --^=4 = -(h + b 2 )r - (B 1 + B 2 ) (3) 

with the real constants aj, Aj, bj, Bj parametrizing the family. As has been re- 
peatedly explained in Refs. [2-5], parameters A 12 of solution (3) are inessential 
for physics in question and we can consider A\ = A 2 = 0. Obviously we have 
YJ^ =1 Ajt = Z)f=i Aj,p = which reflects the fact that for any matrix T from 
SU(3)-Lie algebra we have TrT = 0. Also, as has been repeatedly discussed 
by us earlier (see, e. g., Refs. [2,3]), from the above form it is clear that the 
solution (3) is a configuration describing the electric Coulomb-like colour field 
(components Af' 8 ) and the magnetic colour field linear in r (components A^ s ) 
and we wrote down the solution (3) in the combinations that are just needed 
further to insert into the Dirac equation (4). 

Another part of the family is given by the unique nonperturbative modulo 
square integrable solutions of the Dirac equation in the confining SU(3)-field of 
(3) ^ = (^1, ^2, ^3) with the four- dimensional Dirac spinors ^ representing 
the jth colour component of the meson, so \1/ may describe relative motion 
(relativistic bound states) of two quarks in mesons and the mentioned Dirac 
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equation looks as follows 



where Hamiltonian Hj is 














= H^ = 





H 2 






V o 








(4) 



7 



Ho - i^dr - ij 2 - (do + JtV) - ^7 3 — r— j (fy, + J sin 177V + \ cos 
r V 2 / r sin 1/ V 2 2 



i sini?"^ 



(5) 



with the gauge coupling constant g while /j,q is a mass parameter and one 
should consider it to be the reduced mass which is equal, e. g., for quarkonia, 
to half the current mass of quarks forming a quarkonium. 

Then the unique nonperturbative modulo square integrable solutions of (4) 
are (with Pauli matrix a"i) 



e -iuijt r - 



F j2 (r)a^0^) 



,3 = 1,2,3 



(6) 



with the 2D eigenspinor <&_,- — ^ J of the Euclidean Dirac operator V on 

the unit sphere § 2 , while the coordinate r stands for the distance between 
quarks. The explicit form of <&_,- is discussed in Appendix A. We can call the 
quantity ouj relative energy of jth colour component of meson (while ipj is wave 
function of a stationary state for jth colour component) but we can see that 
if we want to interpret (4) as equation for eigenvalues of the relative motion 
energy, i. e., to rewrite it in the form Hip = cuip with ip = (ipx, ip 2 , ipz) then we 
should put uj = ujj for any j so that Hjipj = ujjipj = cuipj. Under this situation, 
if a meson is composed of quarks q± t2 with different flavours then the energy 
spectrum of the meson will be given by e = m qi + m q2 + u with the current 
quark masses m qk ( rest energies) of the corresponding quarks. On the other 
hand for determination of Uj the following quadratic equation can be obtained 
[1-3] 

[g 2 a 2 +(n j +a j ) 2 ]u 2 -2(\ j -gB j )g 2 a j b j u; j +l(\ j -gB j ) 2 -(n^ 



that yields (at g ^ 0) 



A j g 2 a j b j ± (rij + aj)^(n 2 + 2njaj + A])/ig + g 2 b 2 (n 2 + 2njdj) 

n 2 + 2n i a j + A 2 



(7) 



,j = l,2,3, 
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where g 3 = -(ai + 02), h = -(61 + 6 2 ), B 3 = -(B 1 + 5 2 ), Aj = Aj - 
ctj = \J ftf — g 2 ctj, n-j = 0,1,2,..., while Xj = ±(lj + 1) are the eigenvalues 
of Euclidean Dirac operator Vq on unit sphere with lj = 0, 1,2, .... It should 
be noted that in the papers [1-5] we used the ansatz (6) with the factor e l0Jjt 
instead of e~ l< ^* but then the Dirac equation (4) would look as —id t ^ = H^f 
and in equation (7) the second summand would have plus sign while the first 
summand in numerator of (8) would have minus sign. In the papers [6,7] we 
returned to the conventional form of writing Dirac equation and this slightly 
modified the equations (7)-(8). In the given paper we conform to the same 
prescription as in Refs. [6,7]. 

In line with the above we should have 00 = uji = cj 2 = ^3 in energy spectrum 
e = m qi + m q2 + 00 for any meson (quarkonium) and this at once imposes two 
conditions on parameters dj, bj, Bj when choosing some experimental value for 
e at the given current quark masses m qi ,m q2 . 

The general form of the radial parts of (6) is considered in Appendix B. Within 
the given paper we need only of the radial parts of (6) at rij = (the ground 
state) that are [see (-B.5)] 

F jl = C, I>,r-< ■ ' (l - ^ , P 3 = gb 3 + ft, 

F j2 = K [,()/-< " (l + 9 -^j , Qj = im>- "3 (9) 

with (3j = yj fiQ — ujj + g 2 b 2 while Cj is determined from the normalization con- 
dition f£°(\F jl \ 2 + \F j2 \ 2 )dr = |. Consequently, we shall gain that G L%(R 3 ) 
at any i G R and, as a result, the solutions of (6) may describe relativistic 
bound states (mesons) with the energy (mass) spectrum e. 



3.1 Nonrelativistic limit 



It is useful to specify the nonrelativistic limit (when c — > 00) for spectrum 



(8). For that one should replace g — > g/yfic, a,j — > aj/yhc, bj 
Bj — > Bj/Vhc and, expanding (8) in z = 1/c, we shall get 



bjVhc, 



Uj{nj,lj,\j) = 



±[X Q C 



it 



g 2 a 2 



2h 2 (n j + I A, 



+ 



Xjtnijbj _ g 3 B j a 2 f(n j , Xj) 



M n 3 + \ x 3-\y 



=F Mo- 



h*(n j + \Xj\y 



where f( nj , A,) = AX.n^n 2 + X 2 ) + ^ (raj + 6n 2 X 2 + Xfj. 



z+0(z 2 ), 
(10) 
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As is seen from (10), at c — > oo the contribution of linear magnetic colour 
field (parameters bj, Bj) to spectrum really vanishes and spectrum in essence 
becomes purely nonrelativistic Coulomb one (modulo the rest energy). Also it 
is clear that when rij — > oo, ujj — > ±yVo + 9 2 ^j- 

We may seemingly use (8) with various combinations of signes (±) before sec- 
ond summand in numerators of (8) but, due to (10), it is reasonable to take 
all signs equal to plus which is our choice within the paper. Besides, as is not 
complicated to see, radial parts in nonrelativistic limit have the behaviour of 
form Fj 1 ,Fj 2 ~ r^' +1 , which allows one to call quantum number lj angular 
momentum for jth colour component though angular momentum is not con- 
served in the field (3) [1,3]. So for meson (quarkonium) under consideration 
we should put all lj — 0. 

3.2 Eigenspinors with A = ±1 

Finally it should be noted that spectrum (8) is degenerated owing to degen- 
eracy of eigenvalues for the Euclidean Dirac operator T>q on the unit sphere 
§ 2 . Namely, each eigenvalue of T>q A = ±(/ + 1),/ = 0,1,2..., has multiplicity 
2(/ + 1) so we has 2(1 + 1) eigenspinors orthogonal to each other. Ad referen- 
dum we need eigenspinors corresponding to A = ±1 (/ = 0) so here is their 
explicit form [see (A 16)] 

with the coefficient C = 1/\/2tt (for more details see Appendix A). 

4 Electric form factor, the root-mean-square radius and magnetic 
moment 

As has been mentioned in Section 1, at present little is known about exper- 
imental spectroscopy of the toponium so we should choose a few quantities 
that are the most important from the physical point of view to characterize the 
toponium and then we should evaluate the given quantities within the frame- 
work of our approach. Under the circumstances let us settle on the ground 
state energy of toponium, the root-mean-square radius of it and magnetic mo- 
ment. All three magnitudes are essentially nonperturbative ones and can be 
calculated only by nonperturbative techniques. 
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Within the present paper we shall use relations (8) at rij — = lj so the 
ground state energy of toponium is given by e = 2m t + u with u = u>j(0, 0, Xj) 
for any j = 1, 2, 3 whereas 

g 2 a 1 b 1 a^o g 2 a 2 b 2 a 2 fi g 2 a>zh a 3 fi 

00 = — i ^ nrr = — i ^ TITT = — i ^ t* — r = e ~ 2m t ( 12 ) 

Ai |Ai| A 2 |A 2 | A 3 |A 3 | 

and, as a consequence, the corresponding toponium wave functions of (6) are 
represented by (9) and (11). As the concrete value of e we shall take the one 
of Ref. [14] equal to 347.4 GeV. 



4-1 Choice of t- quark mass and gauge coupling constant 



It is evident for employing the above relations we have to assign some values to 
t-quark mass and gauge coupling constant g. In accordance with Ref. [15] we 
take m t = 173.25 GeV at present. Under the circumstances, the reduced mass 
/la of (5) will take value m t /2. As to gauge coupling constant g = \/4na s , it 
should be noted that recently some attempts have been made to generalize the 
standard formula for a s = a s (Q 2 ) = 127r/[(33 — 2nf) In (Q 2 /A 2 )] (n/ is number 
of quark flavours) holding true at the momentum transfer \fQ 1 — > oo to the 
whole interval < yfQ 1 < oo. We shall employ one such a generalization used 
in Refs. [16]. It looks as follows (x = \J~Cp in GeV) 

l \ 127T f x {x) 

a[X) ~ (33-2n / )l I1 £!±^M 1 6) 

with 

fi(x) = l+(^ {1+X) f 2 ' 2nf) m^-l) 1 + 0.6x 13 j ,f 2 (x) =m 2 (l+2.8x 2 ) 

wherefrom one can conclude that a s — > 7r = 3.1415... when x — > 0, i. e., 
g _> 2tt = 6.2831.... We used (13) at m = 1 GeV, A = 0.234 GeV, n f = 6, 
x = 2mt = 346.50 GeV to obtain g = 1.243528161 necessary for our further 
computations at the mass scale of toponium. 



4-2 Electric form factor 



For each meson (quarkonium) with the wave function ^ = (^j) of (6) we can 
define electromagnetic current = *(/ 3 <g> 7^)^ = (^ t *,* t (/ 3 ® a)^) = 
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(p, J), a = 7°7. Electric form factor f(K) is the Fourier transform of p 
f(K) = [ tfVe-^fx = £ / Vfre^cPx = j^MK) = 

J 7=1 J .7=1 



-iKr 



/e 
(|Fji| 2 + |Fj2| 2 )$j$j~ 2-d 3 x, rf 3 x = r 2 sintfdr<Md<p (14) 

with the momentum transfer X. At rij = = Zj, as is easily seen, for 
any spinor of (11) we have = l/(47r), so the integrand in (14) does 

not depend on tp and we can consider vector K to be directed along z- 
axis. Then Kr = Kr cosd and with the help of (9) and relations (see Ref. 



[17]): f-r^e-^dr = r(a)p-, Re a,p > 0, J^r^e^ (^[^ 



-2 , „2\- a /2 ( sin (a arctan (fT/p)) ^ 



Y(a)(K 2 +p 2 )- a l 2 ) I )J, \\ , Re a > -1, Rep > llmXl, T(a + 
V cos (a arctan (K /p)) / 

1) = aiT (a), Jo e- iKrcos *smi)dtf = 2 sin (Kr)/(Kr), we shall obtain 

» ' (2ft)*** 1 sin [2a 3 - arctan (/r/(2ft))] 

^ M j ^ 6a, ' K(K 2 + 4/3 2 )^ 

wherefrom it is clear that f(K) is a function of i^ 2 , as should be, and we can 
determine the root-mean-square radius of meson (quarkonium) in the form 



< r >- 



* 2a] + 3a j + 1 

\£ w ' (16) 



When calculating (15) also the fact was used that by virtue of the normal- 
ization condition for wave functions we have C 2 [P 2 (1 — gbj/Pj) 2 + Q](l + 

^) 2 ] = (2^ +1 /[3r(2«, + i)]. 

It is clear, we can directly calculate < r > in accordance with the standard 

quantum mechanics rules as < r >= \J J r 2 ^/^d 3 x = \j I r 2 ^^jd 3 x and 

the result will be the same as in (16). So we should not call < r > of (16) the 
charge radius of meson (quarkonium)- it is just the radius of meson (quarko- 
nium) determined by the wave functions of (6) (at rij — = lj) with respect 
to strong interaction, i.e., radius of confinement. Now we should notice the 
expression (15) to depend on 3-vector K. To rewrite it in the form holding 
true for any 4- vector Q, let us remind that according to general considerations 
(see, e.g., Ref. [18]) the relation (15) corresponds to the so-called Breit frame 
where Q 2 = —K 2 [when fixing metric by (1)] so it is not complicated to rewrite 
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(15) for arbitrary Q in the form 



/WJ-jLAW)-!, 6aj ■ ^ w _ g2) „, (17) 

which passes on to (15) in the Breit frame. 



4-3 Magnetic moment 



We can define the volumetric magnetic moment density by m = q(r x J)/2 = 
q[(yJ z — zJ y )i + (zJ x ~ x Jz)i + {xJy — yJ x )k]/2 with the meson charge q and 
J = \l/t(/ 3 (g) Q/)^r. Using (6) we have in the explicit form 

3 <|>t(J). 3 $^<72(Xi$- 
J x = J2( F jl F J2 + F j2 F jl)^T, Jy = J2( F jl F V ~ F j2 F 3l)^—^ "> 

J=l ' J=l ' 

J'=l 

with Pauli matrices Ci,2,3- Magnetic moment of meson (quarkonium) is M = 
md 3 x, where V is volume of meson (quarkonium) (the ball of radius < r >). 
Then at rij = lj = 0, as is seen from (9), (11), Fj\ = F jl ,F* 2 = -F j2 , 

&j(j2Vi®j = for any spinor of (11) which entails J x = J y = 0, i.e., m z = 
while J v m Xy yd 3 x = because of turning the integral over (p to zero, which is 
easy to check. As a result, magnetic moments of mesons (quarkonia) with the 
wave functions of (6) (at lj = 0) are equal to zero, as should be according to 
experimental data [15]. 

Though we can also evaluate magnetic form factor F(Q 2 ) of meson (quarko- 
nium) as a function of Q 2 (see Refs. [6,7]) but the latter will not be used in 
the given paper so we shall not dwell upon it. 



5 An estimate of < r > from leptonic width 



The question now is how to estimate < r > independently to then calculate it 
within framework of our approach. For this aim we shall employ the possible 
width of leptonic decay V — > e + e" which is approximately equal to i\ ~ 13 
keV according to Ref. [13] and V stands for the toponium state analogous 
to J/ip state in charmonium. Under this situation one can use a variant of 
formulas originating from Ref. [19]. Such formulas are often employed in the 
heavy quarkonia physics (see, e. g., Ref. [20]). In their turn they are actually 
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based on the standard expression from the elementary kinetic theory of gases 
(see, e. g., Ref. [21]) for the number v of collisions of a molecule per unit time 



v = y/2 a <v>n, (19) 

where a is an effective cross section for molecules, < v > is a mean molecular 
velocity, n is the concentration of molecules. If replacing v — > r\ we may fit 
(19) to estimate the leptonic width I\ = ri(V — > e + e~) when interpreting 
a as the cross section of creation of e + e~ from the pair it due to electro- 
magnetic interaction, < v > and n as, respectively, a mean quark velocity 
and concentration of quarks (antiquarks) in toponium. To obtain a in the 
explicit form one may take the corresponding formula for the cross section of 
creation of e + e~ from the muon pair (see, e. g., Ref. [18]) and, after 

replacing a em — > Qa em , m M — > m t with electromagnetic coupling constant 
a errt =l/137. 0359895 and muon mass m M , obtain 

> = ^( 1 + ^)^f, (20) 

where electron mass m e = 0.510998918 MeV, the Mandelstam invariant s = 
2m t (m t + e/2) with e from (12), iV is the number of colours and Q = 2/3 
for toponium. To get < v > one may use the standard relativistic relation 
v = \JT{T + 2E )/ (T + E ) with kinetic T and rest energies E for velocity v 
of a point-like particle. Putting T = e/2 — m t , E = m t we shall gain 



< v >= 




At last, obviously, n = 1/V, where the volume of quarkonium V = 4n < r > 3 
/3 with the sought < r >, the latter being yet not related to formula (16). 
The relations (19)-(21) entail the sought independent estimate for < r > 



< r >= ( • j < 22) 

with a of (20). When inserting iV = 3, e = 347.4 GeV, m t = 173.25 GeV, m e = 
0.510998918 MeV, r\ = 13 keV into (22) we shall have < r >^ 0.2162653913 x 
10~ 2 fm. In further considerations we can use this independent estimate of 
< r > while calculating < r > according to (16) which will impose certain 
restrictions on parameters of the confining SU(3)-gluonic field in toponium. 

It should be noted that in the heavy quarkonia physics (see, e. g., Ref. [20]) 
in (19) one often puts n = |-^(0)| 2 , where ^(r) is a wave function of the heavy 
quarkonium stationary state which may be obtained, for example, within 
the framework of potential approach as a solution of the Schrodinger type 
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equation. Following this prescription in our approach with wave functions 

3 

E 



3 

of (6) we should put n = E 1^(0) | 2 with ipj of (6) which would entail 



3 

< r >= [3/(4vr E |^j(0)| 2 )] 1/3 instead of (16). But it is clear that (16) gives 

physically more correct expression for < r > since it employes all values of 
meson wave function rather than the only one at r = (inasmuch as here 



< r >= yj J r 2 ip^ipd?x = J2 I r 2 \vl)j\ 2 d 3 x). So we shall use just (16) in what 
follows. 



6 Estimates for parameters of SU(3)-gluonic field in the ground 
state of toponium rj t 



Now we are able to estimate parameters a,j, bj, Bj of the confining SU(3)-field 
(3) for the toponium ground state r\ t within framework of two approaches - 
relativistic and nonrelativistic ones. 



6.1 Relativistic approach 



Under this approach we should consider (12) and (16) the system of equations 
which should be solved compatibly if taking e = 347.4 GeV, m t = 173.25 
GeV and < r >m 0.2162653913 x 10~ 2 fm in accordance with the independent 
estimate of Section 5. While computing for distinctness we take all eigenvalues 
Xj of the Euclidean Dirac operator V on the unit two-sphere § 2 equal to (- 
1). The results of numerical compatible solving of equations (12), (16) are 
adduced in Tables 1-2. 

Table 1 

Gauge coupling constant, reduced mass hq and parameters of the confining SU(3)- 
gluonic field in the toponium ground state r] t : relativistic approach 



9 


Mo (GeV) 


a 1 


02 


6i (GeV) 


b 2 (GeV) 


Bi 


B 2 


1.24353 


86.6250 


0.361253 


0.339442 


48.9402 


76.7974 


-0.360 


-0.295 



Table 2 

Theoretical ground state energy of toponium and its radius: relativistic approach 



Theoret. e (GeV) 


Theoret. < r > (fm) 


e = 2m t + Wj (0,0,-1) = 347.400 


0.213915 x 10-2 
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6.2 Nonrelativistic approach 



If estimating the quark velocity in toponium by the relation (21) then at 
e = 347.4 GeV, m t = 173.25 GeV we obtain < v >^ 0.071935 that points 
out a nonrelativistic approach to be applicable. It should be noted, however, 
this nonrelativistic approach should be consistent with the above relativistic 
one. We cannot, therefore, follow the standard strategy of the heavy quarkonia 
physics which exploits the so-called potential approach (see, e. g., Refs. [8,20]). 
The essence of the latter is in that the interaction between quarks is modelled 
on a nonrelativistic confining potential in the form V(R) = a/R + kR + cq 
with some real constants a, k, cq and the distance between quarks R. However, 
parameters of such a potential, i.e. quantities a, k, c are not related with QCD- 
Lagrangian in any way and we cannot speak about V(R) as describing some 
gluon configuration between quarks. It would be possible if the mentioned 
potential were a solution of Yang-Mills equations directly derived from QCD- 
Lagrangian since, from the QCD-point of view, any gluonic field should be 
a solution of Yang-Mills equations (as well as any electromagnetic field is by 
definition always a solution of Maxwell equations). 

In reality, as was shown in Refs. [2,3] (see also Appendix C), potential of form 
a/R + kR + Co cannot be a solution of the Yang-Mills equations if simultane- 
ously a 7^ 0, k 0. Therefore, it is impossible to obtain compatible solutions 
of the Yang-Mills-Dirac (Pauli, Schrodinger) system when inserting poten- 
tial of the mentioned form into Dirac (Pauli, Schrodinger) equation. So, we 
draw the conclusion (mentioned as far back as in Refs. [4]) that the potential 
approach seems to be inconsistent: it is not based on compatible nonperturba- 
tive solutions for the Dirac- Yang-Mills system derived from QCD-Lagrangian 
in contrast to our confinement mechanism. Actually potential approach for 
heavy quarkonia has been historically modeled on positronium theory. In the 
latter case, however, one uses the unique modulo square integrable solutions of 
Dirac (Schrodinger) equation in the Coulomb field [condensate of huge number 
of (virtual) photons], i. e., one employs the unique compatible nonperturba- 
tive solutions of the Maxwell-Dirac (Schrodinger) system directly derived from 
QED-Lagrangian to describe positronium (or hydrogen atom) spectrum. 

On the other hand, as was mentioned in Section 1, our confinement mechanism 
is based on the unique family of compatible nonperturbative solutions for the 
Dirac- Yang-Mills system derived from QCD-Lagrangian and just magnetic 
colour field of solution (3) is responsible for linear confinement. But, as we 
have seen in Section 2 [see relation (10)], if directly taking the nonrelativistic 
limit c — > oo then the contribution of linear magnetic colour field (param- 
eters bj,Bj) to spectrum really vanishes and spectrum in essence becomes 
purely nonrelativistic Coulomb one (modulo the rest energy). Consequently, 
as we emphasized as far back as in Refs. [4], the confinement mechanism un- 
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der discussion is essentially connected with relativistic effects conditioned by 
availability of the mentioned magnetic colour field between any two quarks. 
Under the circumstances the only reasonable way of constructing a nonrela- 
tivistic approach within our confinement scheme is the power series expansion 
of the physical magnitudes of interest in z — 1/c with retaining necessary 
number of terms. It is clear, we then shall obtain the consistent transition 
from relativistic regime to nonrelativistic one. In their turn, the mentioned 
magnitudes may be computed within the relativistic framework with the help 
of wave functions (6) and then the necessary expansions should be fulfilled. 

Following the just formulated receipt for description of toponium in a nonrel- 
ativistic manner we should use (10) (at h — c — 1) with subtracting the rest 
energy /i c 2 to replace the relations (12) (at rij = 0, Xj = —1) by 

w = -fi g 2 a 2 1 /2 - g 2 a 1 b 1 + /j l0 g 3 B 1 al = -/j, g 2 a 2 2 /2 - g 2 a 2 b 2 + ^ g z B 2 a 2 2 = 

-fi g 2 a 2 3 /2 - g 2 a 3 b 3 + fi g 3 B 3 a 2 3 = e - 2m t (23) 
with a 3 = -(ai + a 2 ), h = -(bi + b 2 ), B 3 = -(B 1 + B 2 ). 

At the same time we can use (16) to compute < r > with replacing the quanti- 
ties (Xj = yj (Xj — gBj) 2 — g 2 a 2 , (3j = — uj 2 + g 2 b 2 by their nonrelativistic 
expressions 

<*i = W-^9B j z + 0{#), (24) 

Pi = ^ + ^ (g 2 b 2 - ^oA 4 ) z2 + ■ (25) 

Before expanding f3j we made replacement Uj — > ujj — fi c 2 in the formula 
j3j = ^J/jLq — uj 2 + g 2 b 2 , i. e. we subtracted the rest energy from uj as is required 
in nonrelativistic limit (see, e. g., Ref. [18]). After it we should compatibly 
solve equations (23) and (16) with aj,/3j of (24)-(25) (with % — c — 1) at 
Xj = -1, e = 347.4 GeV, m t = 173.25 GeV and < r >^ 0.2162653913 x 10~ 2 
fm in accordance with the independent estimate of Section 5. When solving 
we should impose the conditions otj > —1/2, (3j > to have ^ e L 2 (R 3 ) at 
any t G 1 for ^ of (6). 

The results of numerical computation are adduced in Tables 3-4. 
Table 3 

Gauge coupling constant, reduced mass [1q and parameters of the confining SU(3)- 
gluonic field in the toponium ground state 7ft : nonrelativistic approach 



9 


/i (GeV) 


a x 


02 


h (GeV) 


b 2 (GeV) 


Bi 


B 2 


1.24353 


86.6250 


0.200 


1.32600 


-291.248 


143.349 


-0.780 


1.40873 
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Table 4 

Theoretical ground state energy of toponium and its radius: nonrelativistic approach 



Theoret. e (GeV) 


Theoret. < r > (fm) 


e = 2m t + u)j(0, 0, -1) = 347.400 


0.216264 x 10~ 2 



7 Consistency with possible width of two-photon decay rj t — > 27 



Let us consider whether the estimates of previous section are consistent with 
possible width T 2 of two-photon decay r\ t — > 27 which might be one of the 
main signals for detecting toponium under certain conditions (see Section 2). 
To estimate T 2 we can use an analogy of toponium with charmonium where 
width T(i] c — > 27) ~ 7.2 keV for the charmonium ground state 7] c whereas 
leptonic width T(J/ip — > e + e~) ps 5.55 keV for state J/^ [15]. Accordingly, 
ratio V(r] c — > 2^)/r(J/t/j — > e + e~) ps 1.297 and if taking the same ratio for 
r 2 /Ti with Ti ps 13 keV of Section 5 for toponium state V analogous to J/ip 
then we shall obtain T 2 ~ 16.865 keV. 

On the other hand, actually kinematic analysis based on Lorentz- and gauge 
invariances gives rise to the following expression for width T of the electro- 
magnetic decay P — > 27 (where P stands for any meson from ir°, 77, r]', see, 
e.g., Refs. [22]) 

r = ^ira 2 em g 2 Pi y (26) 

with electromagnetic coupling constant a em = 1/137.0359895 and P-meson mass 
H while the information about strong interaction of quarks in P-meson is en- 
coded in a decay constant gp 77 . Making replacement gp 77 = /p//j we can 
reduce (26) to the form 

r = Kal^pp ^ 

Now it should be noted that the only invariant which fp might depend on is 
Q 2 = fi 2 , i. e. we should find such a function J-'iQ 2 ) for that T{Q 2 = fi 2 ) = fp. 
It is obvious from physical point of view that T should be connected with 
electromagnetic properties of P-meson. As we have seen above in Section 4, 
there are at least two suitable functions for this aim - electric and magnetic 
form factors. But there exist no experimental consequences related to magnetic 
form factor at present whereas electric one to some extent determines, e. g., 
an effective size of meson (quarkonium) in the form < r > of (16). It is 
reasonable, therefore, to take T{Q 2 = /i 2 ) = Af(Q 2 = fi 2 ) with some constant 
A and electric form factor / of (17) for the sought relation. Under the situation 
we obtain additional equation imposed on parameters of the confining SU(3)- 
gluonic field in P-meson which has been used in Refs. [6,7] for to estimate 
the mentioned parameters in 7r°- and ^-mesons. Inasmuch as relation (27) 
is in essence nonperturbative since decay constant fp cannot be computed 
by perturbative techniques, we may extend (27) over heavy quarkonia states 
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similar to 7T°, 77, r]', in particular, over r) t . As a result, using (17) we come from 

(27) to relation 

r = T 2 = ^ (a f -J- ■ Sin { p 2 w 16,865 keV (28) 

4 ^ fl (in,,-, (l-x 2 )^ J 1 J 

with Xj = fj,/(2/3j) and fi = e = 347.4 GeV. Under the circumstances we 
can employ the results of Tables 1 and 3 and compute the left-hand side of 

(28) in relativistic and nonrelativistic regimes respectively which entails the 
corresponding values A 0.0340 and A 0.156. Consequently, we draw the 
conclusion that parameters of the confining SU(3)-gluonic field in toponium 
from Tables 1 and 3 might be consistent with r 2 in both regimes. 



8 Estimates of gluon concentrations, electric and magnetic colour 
field strengths 



Now let us remind that, according to Refs. [3,5], one can confront the field (3) 
with Too-component (volumetric energy density of the SU(3)-gluonic field) of 
the energy-momentum tensor (2) so that 

E 2 + H 2 1 ( aj + a^ + al bj + hh + bj \ = A _B_ 
00 _ " 2 2 ^ r 4 + r 2 sin 2 tf ) ~ r 4 + r 2 sin 2 ■& 

(29) 

with electric E and magnetic H colour field strengths and with real A > 0, B > 
0. One can also introduce magnetic colour induction B = (An x 10~ 7 H/m) H, 
where H in A/m. 

To estimate the gluon concentrations we can employ (29) and, taking the 
quantity ou — F± — Fi(V — > e + e~) w 13 keV of Section 5 for the characteristic 
frequency of gluons, we obtain the sought characteristic concentration n in 
the form 

Too / Qn ^ 
n = — (30) 

1 1 

so we can rewrite (29) in the form T o = Tqq 111 + Tqq 1 conforming to the con- 
tributions from the Coulomb and linear parts of the solution (3). This entails 
the corresponding split of n from (30) as n = n cou \ + n\ in . 

The parameters of Tables 1 and 3 were employed when computing and for 
simplicity we put sin-# = 1 in (29), whereas the Bohr radius a = 0.529177249- 
10 5 fm [15]. 

Table 5 contains the numerical results for n cou i, n\ [n , n, E, H, B for the quarko- 
nium under discussion in relativistic approach while Table 6 is obtained in 
nonrelativistic one. 
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Table 5 



Gluon concentrations, electric and magnetic colour field strengths in toponium: 
relativistic approach 



'IX ■ 


r =< r > = 0.213915 x 10~ 2 


fm 




















r (fm) 


T} i ( in 3 ^ 






n (m~ 3 ) 


E (V /m) 




H (A/m) 




B (T) 






0.1r 


0.880566 X 10 66 




0.131833 X 1() 64 


0.881884 X 10 66 


0.223572 


X 10 30 


0.116367 x 


1() 27 


0.146231 


X 


10 21 


T-0 


0.880566 X 10 62 




0.131833 X 1() 62 


0.101240 X 10 63 


0.223572 


X 10 2S 


0.116367 x 


1Q 26 


0.146231 


X 


10 20 


10r 


0.880566 X 10 58 




0.131833 X 10 60 


0.140639 X 10 60 


0.223572 


X 10 26 


0.116367 x 


10 25 


0.146231 


X 


10 19 


1.0 


0.184386 X 10 52 




0.603264 x 10 56 


0.603282 x 10 56 


0.102305 


X 10 23 


0.248927 x 


10 23 


0.312811 


X 


10 17 


a 


0.235138 X 10 33 




0.215429 x 10 47 


0.215429 x 10 47 


0.365339 


X 10 13 


0.470404 x 


1() 1S 


0.591127 x 


10 12 


Table 6 
























Gluon concentrations, electric 


and 


magnetic colour field strengths 


in toponium: 










nonrelativistic approach 






















Vf- 


r =< r > = 0.216264 x 10~ 2 


fm 




















r (fm) 


"coul (m _3 ) 




n Mn (m 3 ) 


n (m~ 3 ) 


E (V/m) 




H (A/m) 




B (T) 






O.lro 


0.402970 X 10 67 




0.680959 X 10 64 


0.403651 x 10 67 


0.478269 


X 10 30 


0.264472 x 


1() 27 


0.332345 


X 


10 21 


ro 


0.402970 X 10 63 




0.680960 X 10 62 


0.471066 x 10 63 


0.478269 


X 10 2S 


0.264472 x 


1() 26 


0.332345 


X 


10 20 


l()r 


0.402970 X 10 59 




0.680960 X 10 60 


0.721257 x 10 60 


0.478269 


X 10 26 


0.264472 x 


10 25 


0.332345 


X 


10 19 


1.0 


0.881475 X 10 52 




0.318486 X 10 57 


0.318494 x 10 57 


0.223687 


X 10 23 


0.571957 x 


10 23 


0.718742 




10 17 


ao 


0.112410 X 10 34 




0.113733 X 10 4S 


0.113733 x 10 48 


0.798800 


X 10 13 


0.108084 x 


10 19 


0.135823 


X 


10 13 



9 Discussion and concluding remarks 



9. 1 Discussion 

We can see that main reasonable characteristics of r\ t such as energy of ground 
state and the root-mean-square radius (in essence, radius of confinement) and 
also width of possible decay rjt — > 27 may be consistent with appropriate pa- 
rameters of the confining SU(3)-gluonic field between quarks in toponium in 
both relativistic and nonrelativistic approaches within the framework of our 
confinement mechanism. In other words, we can obtain a description of to- 
ponium directly appealing to quark and gluonic degrees of freedom as should 
be from the first principles of QCD. As is seen from Tables 5 and 6, at both 
relativistic and nonrelativistic description the gluon concentrations are huge 
at the characteristic scales of the toponium ground state f] t and the corre- 
sponding fields (electric and magnetic colour ones) can be considered to be 
the classical ones with enormous strengths. The part n cou i of gluon concen- 
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tration n connected with the Coulomb electric colour field is decreasing faster 
than nun, the part of n related to the linear magnetic colour field, and at large 
distances n\ in becomes dominant. It should be emphasized that in fact the 
gluon concentrations are much greater than the estimates given in Tables 5 
and 6 because the latter are the estimates for sufficiently big possible gluon 
frequencies, i.e. for sufficiently big possible gluon impulses (under the concrete 
situation of toponium ground state rjt). As was mentioned in Section 1, the 
overwhelming majority of gluons between quarks should be soft, i. e., with 
frequencies << 13 keV so the corresponding concentrations >> the ones in 
Tables 5 and 6. The given picture is in concordance with the one obtained in 
Refs. [4-7]. As a result, the confinement mechanism developed in Refs. [1-3] 
is also confirmed by the considerations of the present paper. 

By the way, the estimate of < r >m 0.216 x 10~ 2 fm obtained in present pa- 
per allows one to make a suggestion about why observation of toponium finds 
difficulty at colliders. Let us use an analogy with classical electrodynamics 
where, as is well known (see e. g. Ref. [23]), the notion of classical electromag- 
netic field (a photon condensate) generated by a charged particle is applicable 
only at distances >> the Compton wavelength A c = 1/m for the given point- 
like particle with mass m. Passing on to QCD, gluons and quarkonia and 
replacing electromagnetic field by colour one in the case of t-quark with mass 
m t = 173.25 GeV we obtain A c xs 0.11 x 10~ 2 fm. If comparing the above < r > 
to characteristic radius of weak interaction r wea k ~ 1/irty/ m 0.2450 x 10~ 2 fm 
(mw ~ 80.403 Gev is mass of W^-boson) and to the just obtained A c then we 
have inequality A c < r wea k ~ < r > so one may draw the conclusion that in 
pair it when creating at colliders the most probable distance between quarks 
is < r wea k so that quarks are more inclined to weak interaction rather than to 
strong one. In other words, they have not time in order to form a classical con- 
fining SU(3)-gluonic field due to strong interaction and to constitute a bound 
state in virtue of it. But we cannot completely exclude events where distance 
between t-quarks when creating at colliders would be much greater than A c 
which might entail formation of r) t , for example, and, consequently, a signal 
for detecting toponium, e.g., in the form of decay r) t — > 27, as mentioned in 
Section 2. 

It should be noted, however, that our results are of a preliminary character 
not only because of that the experimental spectroscopy of toponium is in its 
infancy but also in virtue of that the current quark masses (as well as the gauge 
coupling constant g) used in computation are known only within the certain 
limits and we can expect similar limits for the magnitudes discussed in the 
paper so it is necessary further specification of the parameters for the confining 
SU(3)-gluonic field in toponium which can be obtained when experimental 
situation for toponium becomes more satisfactory. We hope to then continue 
analysing the toponium physics. 
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9.2 Concluding remarks 



Finally we should note the following. As has been shown in the paper, our 
approach allows one to conduct both relativistic and nonrelativistic description 
and both the cases are consistent with each other. Only experiments can, 
however, determine what physical picture (relativistic or nonrelativistic one) 
for quarkonia is really realized and enough for their complete description. 
Our approach works in either case since it is based on the unique family of 
compatible nonperturbative solutions for the Dirac- Yang-Mills system derived 
from QCD-Lagrangian and, as a result, the approach is itself nonperturbative, 
relativistic from the outset, admits self-consistent nonrelativistic limit and may 
be employed for any meson (quarkonium). 



Appendix A 



We here represent some results about eigenspinors of the Euclidean Dirac 
operator on two-sphere S 2 employed in the main part of the paper. 

When separating variables in the Dirac equation (4) there naturally arises 
the Euclidean Dirac operator T>q on the unit two-dimensional sphere § 2 and 
we should know its eigenvalues with the corresponding eigenspinors. Such a 
problem also arises in the black hole theory while describing the so-called 
twisted spinors on Schwarzschild and Reissner-Nordstrom black holes and it 
was analysed in Refs. [3,24], so we can use the results obtained therein for our 
aims. Let us adduce the necessary relations. 

The eigenvalue equation for corresponding spinors $ may look as follows 

£> $ = A$. (Al) 



As was discussed in Refs. [24], the natural form of V in local coordinates f 
on the unit sphere § 2 looks as 



V Q = -io x 



1 

'sin?? 



io- 2 d$ + ia-i-r— - ( d v - -0203 cosfi 



1 cot i? 

o\o 2 d^ + - — -0x0$^ — 0-102 

sin v 2 



with the ordinary Pauli matrices 



(A2) 



o-i = 



1 

1 



,02 = 



-i 

1 



,03 = 



1 

-1 
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so that o{Dq = —T>q<j\. 



The equation (Al) was explored in Refs. [24]. Spectrum of D consists of the 
numbers A = ±(7+1) with multiplicity 2(/ + l) of each one, where / = 0, 1, 2, .... 
Let us introduce the number m such that — / < m < Z+l and the corresponding 
number m' = m — 1/2 so \m'\ < I + 1/2. Then the conforming eigenspinors of 
operator V are 

*= (£)=**= ? ( p k in % ?" /2 ) e ~' m " f (A3) 

with the coefficient C = and k = 1 + 1/2. These spinors form an or- 

thonormal basis in L^S 2 ) and are subject to the normalization condition 

y$ t $^=y yd^xi 2 + i$ 2 i 2 )smw^ = i . (aa) 

S 2 

Further, owing to the relation aiV Q = —T>qO\ we, obviously, have 

(7i$ TA = $ ±A . (A5) 



As to functions P^, n , (cos = P^, n ,(cos^) then they can be chosen by mis- 
cellaneous ways, for instance, as follows (see, e. g., Ref. [25]) 



P*, n ,(cos#) 



m'+n' 
2 



(k - m')\(k - n')\ (\ + costf' 
\j {k + m')\(k + n')\ ^1 -cos^ 

(A; + j)H 2i /l-costfV 



x 



x y 

. i , (A;- j)!0'-"i') ! 0'-«') ! 

j=max(m ,n ) v y VJ 7 



with the orthogonality relation at m', n' fixed 



2 

| P^(cos^)P^ n ,(cos^)sin«= 2k+~l 5kk ' 



(A6) 



(A7) 



It should be noted that square of V is 

COS $ 

£> 2 = -A s 2/ 2 + cr 2 or 3 . 2 Q d v + 

sin v 



+ 



4 sin 2 d 4 ' 



while Laplacian on the unit sphere is 



(A8) 



sin-# 



d# sin W# + — -n-;dl = dl + cot ^ + -r-o-j^ 
sin 17 ' , . - < 



sin 2 ?? 



(A9) 



so the relation (A8) is a particular case of the so-called Weitzenbock-Lichnerowicz 
formulas (see Refs. [26]). Then from (Al) it follows £> 2 $ = A 2 $ and, when 
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using the ansatz $ = P(^)e" imV = (p 1 ) e"^, P 1;2 = Pi, 2 (tf), the equation 



p2<£ _ turns j nto 



m' 2 + 7 m' cos $ 



-<9 2 - cot W e + , 2 ~ + "".7V ^ 1 P = 

A 2 -1)P, (A10) 
wherefrom all the above results concerning spectrum of V can be derived [24]. 

When calculating the functions P^ l , n ,{cos , &) directly, to our mind, it is the 
most convenient to use the integral expression [25] 



P*, n ,(cos0) = — « 



(k - m')\(k + m')\ 



2n 



2tt\ (k - n')\(k + ri) 



0! J 



iip/2 $ 

3»V/ <■ one — 



cos — + ie 



-^ 2 sin- 



k—n' 



X 



^/ 2 sin? + e-^ 2 cos^ 



k+n' 



im (p 



dip 



2 2 / 
and the symmetry relations (z = cos $) 

Pm'n'( Z ) = Pn'm'( Z )' ^m> ,-n'( Z ) = P-m\n'( Z )i Pm'n'( Z ) = P-m',-n'( Z ) > 

r>k ( \ ■2k-2m'~2n / r>k 

r m'n'\~ Z ) — 1 

In particular 



(All) 



Pm'-n'( Z ) ■ 



(A12) 



i k \/(2kV 

P k kki z ) = cos 2fc (0/2), PLki z ) = ? k sm 2fe (0/2), PU Z ) = A^T 1 sm fc * , 



2 fe A;! 



■k—n 



(2k)\ 



— sin fe " n ' (0/2) cos k+n ' (0/2) . (A 13) 



\| (ife-nOKifc + n')! 



If A = ±(Z + 1) = ±1 then I = and from (A3) it follows that fc = I + 1/2 = 

1 /2 

1/2, |m'| < 1/2 and we need the functions P ±1 / 2 ±i/ 2 that are easily evaluated 
with the help of (All)-(A13) so the eigenspinors for A = — 1 are 



<I> =-( C ° S i + tSin l)e^,<S>=C 



2 V cos | — i sm | 



cos | + « sin £ 



2 ) g-^/2 

2 V — cos | + « sin | ' 



(A14) 



while for A = 1 the conforming spinors are 



$ = — 



C / cos I — « sin % \ ,-,„/ 9 ^ C / — cos f + « sin f 



2 V cos | + ? sm | 



2 V cosf + zsmf ' 
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with the coefficient C = yl/(2n). 

It is clear that (A14)-(A15) can be rewritten in the form 

A = -1 : * = % { e %\ e^ 2 , or $ = % ( e % ) e~*"\ 

A = 1: * = l( e e'?) e "" 2 ' " * = 1 ( "Ji 1 ) e "*" /2 ' (A16 > 
so the relation (A5) is easily verified at A = ±1. 



Appendix B 



We here adduce the explicit form for the radial parts of meson wave functions 
from (6). At rij — they are given by 



Y, 



while at rij > by 

F n = CjPjr^e-^ 



F j2 = X '.(J, /■'•■< i;r 



^7 nj 7 ^ 



(P.2) 



with the Laguerre polynomials L^(rj), rj = 2(3jr, (3j = — a;| + g 2 & 2 at 

j = 1, 2, 3 with 6 3 = -(61 + 6 2 ), Pj = + <5j = Mo - ^j, = PjQj®j + 
(P/ - Q]) gaj /2, = Pj-Qj-Aj + (P/ + Q])gaj/2 with a 3 = -(ai + a 2 ), Aj = 

Aj — gBj with P 3 = —(Pi + B 2 ), aj = 



/A] - g 2 a], while A^ = ±(/ j + 1) 
are the eigenvalues of Euclidean Dirac operator T> on unit two-sphere with 
lj = 0,1,2,... (see Appendix A) and quantum numbers rij = 0,1,2,... are 
defined by the relations 

gbjZj - PjYj 



nj = 



PjPjQj 

which entails the quadratic equation (7) and spectrum (8). Further, Cj of 
(P.1)-(P.2) should be determined from the normalization condition 



00 

/(l^-i| 2 + l^ 2 | 2 )^ = -. 



(5.4) 
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As a consequence, we shall gain that in (6) tyj e L^K 3 ) at any i 6 1 and, 
accordingly, ^ = (vE^, \l/ 2 , ^3) m &y describe relativistic bound states in the 
field (3) with the energy spectrum (8). As is clear from (B.3) at rij = we 
have gbj/Pj = Yj/Zj so the radial parts of (B.l) can be rewritten as 



/•}, = CjIY"' 11 ( l - !J J-)- F * = i( ¥hr ni <- 3jr (l • ^) • (5-5) 
More details can be found in Refs. [1,3]. 



Appendix C 

The facts adduced here have been obtained in Refs. [1-3] and we concisely 
give them here only for completeness of discussion in Section 6. 

The Dirac- Yang-Mills system derived from QCD-Lagrangian according to the 
standard prescription of Lagrange approach is 

X^ = /i #, (C.l) 

d*F = g{*F A A- AA*F) + gj (C.2) 

with a gauge coupling constant g, Dirac operator D, F = dA + gA A A and 
the Cartan's wedge (external) product A, whereas * means the Hodge star 
operator conforming to a Minkowski metric, for instance, in the form of (1), 
while the source J (a non-Abelian SU(3)-current) is 

j = j;\ a * (dxn = *j = *(j;\ a dxn = *{j a \a) , (C3) 

where currents 

j a = j;dx» = y(i 3 ®^)\ a ydx», 

so summing over a — 1, ...,8 is implied in (C.3). Besides we have div(j a ) = 
div(j) = if \1/ obeys Dirac equation C.l [1,3], where the divergence of the 
Lie algebra valued 1-form A = A^dx 11 = A a ^\ a dx^ is defined by the relation 
(see, e. g. Refs. [26]) 

div(A) = -Ld,(V5g^A u ) . 

Definitions of the operators * and d (external differentiation) can be found in 
Refs. [26] while explicit form of Dirac operator T> of (C.l) depends on choice of 
local coordinates on Minkowski spacetime and for case of coordinates t, r, 1?, <p 
Dirac equation (C.l) can be rewritten in form (4)-(5) of Section 3 and if we 
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require its modulo square integrable solutions to consist from the components 
of form tyj ~ r aj e~* jr with some otj > 0, f3j > then it will entail all the 
components of the current J to be modulo << 1 at each point of Minkowski 
space (perhaps except for a small neighbourhood of point r = 0). The latter 
allow us to put J ~ and we come to the problem of finding the confining 
solutions for the Yang-Mills equations of (C.2) with J = whose unique 
nontrivial form is given by (3) of Section 3 while the unique corresponding 
modulo square integrable solutions of C.l are given by (6) of Section 3 (for 
more details see Refs. [1-3] and Appendix B). 

As has been mentioned in Subsection 6.2, in meson spectroscopy one often uses 
the nonrelativistic confining potentials. Those confining potentials between 
quarks are usually modelled in the form a/r + br with some constants a and b. 
It is clear, however, that from the QCD point of view the interaction between 
quarks should be described by the whole SU(3)-field = A^X a , genuinely 
relativistic object, the nonrelativistic potential being only some component 
of A$ surviving in the nonrelativistic limit when the light velocity c — > oo. 
Let us explore whether such potentials may be the solutions of the Maxwell or 
SU(3)-Yang-Mills equations. Though this can be easily derived from the above 
result about uniqueness of solution (3) of Section 3, let us consider the given 
situation directly in view of its physical importance. We shall use the Hodge 
star operator action on the basis differential 1- and 2-forms on Minkowski 
spacetime with local coordinates t, r, i?, ip in the form 

*dt = r 2 sin $dr A d*& A d(p, *dr = r 2 sin -&dt A d$ A d(p, 

*d'd = — r sin ddt Adr A dip, *dip = rdt A dr A dti, 
*(dt A dr) = -r 2 sin -ddfi A dip , *(dt A d§) = sin §dr A dip , 

*(dt A dip) = ^—:dr A dd , *(dr A d-d) = sin fldt A dip , 

sin 17 

*(drAd<p) = —dtAdi). (CA) 

sin 17 



Maxwell Equations 

In the case of Maxwell equations (i.e., Yang-Mills equations for the case of 
U(l)-group looking as d* F = at J = 0) the ansatz A = A t dt = (a/r + br)dt 
yields F = dA = (a/r 2 — b)dt A dr. Then with the help of (CA) we have 
*F = sin $(br 2 — a)d'd A dip and the relation d* F = 2br sin $dr A dd A dip = 
entails 6 = 0. 
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SU (3)- Yang- Mills Equations 



We use the ansatz 

A = A a t \ a dt = (A'/r + B'r)dt (C.5) 

with some constant matrices A' = a a \ a , B' = f3 a \ a . Then A A A = 0, F = 
dA + gA A A = dA = {A'/r 2 - B')dt A dr. Again with the help of (C.4) 
we have *F = sm-d(B'r 2 - A') d-d A dip, d * F = 2B'rsm-ddr A d-d A dip, 
*F AA — AA*F = —2[A', B']rdt Adfi Adtp. Under the circumstances the Yang- 
Mills equations (C.2) (at J = 0) are tantamount to the conditions d * F — 0, 
*FAA-AA*F = 0. The former entails 5' = 0, then the latter is fulfilled at 
any A' and we can see that the Coulomb-like field A = (A'/r)dt is a solution 
of the Yang-Mills equations (C.2) (at J = 0) with arbitrary constant matrix 
A'. In principle the ansatz (C.5) might be a solution of (C.2) with the source 
of the form 

J = 2B'r sin Mr A d-d A dip + 2g[A' , B'\rdt A d-d A dip = 

*j = *(j;\ a dxn = * (—dt + g 2 ^ldr\ , (C.6) 
p y r r sin 17 ) 

but div(j) 7^ and this is not consistent with the only source (C.3) derived 
from the QCD-Lagrangian. We can avoid this difficulty putting matrices A', B' 
are not equal to zero simultaneously and both matrices belong to Cartan 
subalgebra of SU(3)-Lie algebra, i.e. commutator [A', B'\ = 0. Then B' = 
+ Ps^8 and for consistency with the only admissible source of (C.3) we 
should require source of (C.3) to be equal to one of (C.6) which entails 

wherefrom one can conclude that 

— — 23 3 

g*(I 3 ® 7 <)A a ^ = 0, a ^ 3, 8 , gV(I 3 ® lt )\ 3 V = — , 

r 

28 8 

g*(h ® lt)\ 8 ^ = -y , g*(h ® 7^)A a ^ = 0, a = 1, 8, n ^ t , (C.7) 

which can obviously be satisfied only at (3 3 ~ j3 8 ~ ^ — > at each point 
of Minkowski spacetime, i. e., really matrix B' — again. All the above can 
easily be generalized to any group SU(iV) with TV > 1 [2,3]. 

As a result, the potentials employed in nonrelativistic approaches do not obey 
the Maxwell or Yang-Mills equations. The latter ones are essentially relativistic 
and, as we can see from (3) of Section 3, the components linear in r of the 
unique exact solution A^ are different from A t and related with magnetic 
(colour) field vanishing in the nonrelativistic limit. 
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Remark about search for nonrelativistic confining potentials 

The above results make us cast a new glance at search of many years for non- 
relativistic potentials modelling the confinement. Many efforts were devoted 
to the latter topic, for example, within the framework of lattice gauge theories 
or potential approach (see, e.g., Ref. [20] and references therein). It should 
be noted, however, that almost in all literature on this direction one does not 
bring up the question: whether such potentials could (or should) satisfy the 
Yang-Mills equations? As is clear from the above the answer is negative. That 
is why the mentioned approaches seem to be inconsistent: such potentials can- 
not describe any gluonic configuration between quarks since any gluonic field 
should be a solution of Yang-Mills equations (as well as any electromagnetic 
field is by definition always a solution of Maxwell equations). 
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